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Abstract. In this presentation we review the link between the statistics of intensity fluctuations 
in spectral line data cubes with underlying statistical properties of turbulence in the interstellar 
medium. Both the formalism of Velocity Channel Analysis for optically thin lines and its extension 
to the lines with self-absorption is described. We demonstrate that by observing optically thin lines 
1 from cold gas in sufficiently narrow (thin) velocity channels one may recover the scaling of the 

stochastic velocities from turbulent cascade, in particular, Kolmogorov velocities give K~ 21 con- 
tribution to the intensity power spectrum. Synthetically increasing the channel thickness separates 
out the underlying density inhomogeneities of the gas. Effects of self absorption, on the other hand, 
\ retain the velocity signature even for integrated lines. As a result, intensity fluctuations tend to show 

C^ - ■ universal but featureless scaling of the power « K~ 3 over the range of scales. 

in 
o 

)Q INTRODUCTION 

O 

There is little doubt that interstellar medium is turbulent (see reviews Jli 0,0]). Turbu- 
lence proved to be ubiquitous in molecular clouds fl], diffuse ionized o] and neutral [6] 
media. Recent years were marked by a substantial progress in theoretical and numerical 
description of both incompressible and compressible Magneto-Hydrodynamical (MHD) 
turbulent cascade O, IM la, [T^, |Tl|, [l^]. MHD turbulence controls many essential astro- 
physical processes including star formation, transport and acceleration of cosmic rays, 
of heat and mass (see H13lll4ll5lll6l.ll0ll ). 

The progress in testing and advancing our theoretical understanding of these astro- 
physical processes rely on recovering the actual parameters of turbulent cascade from 
observational data. 

For this, the statistical approach is useful (see IntE])- In fact, the turbulent medium 
is, generally, a non-uniform distribution of matter with stochastic, random, density 
p(x), that moves according to stochastic velocity field u(x). Recovering the scaling 
of statistical descriptors for stochastic density and velocity will allow us to get insight 
into mechanisms, responsible for turbulent driving. Among the two, the most direct 
information about the turbulence is contained in the statistics of the velocity field, while 
the statistics of density fluctuations usually provides complimentary indirect way of 
testing turbulence. 

A wealth of observational information is contained in spectral or channel maps of 
the emission lines produced by gas tracing the turbulent cascade. Both inhomogeneous 
distribution of the emitting matter and its motions result in small scale fluctuations 
of the emission intensity observed at a given velocity. The problem is to relate the 



two dimensional statistics of intensity maps to underlying three dimensional properties 
of the turbulence and disentangle the effects of random motions from the effects of 
density fluctuations . This problem is addressed in the velocity channel analysis (VCA) 
lfl9l I20I1 which we review here. VCA provides theoretical framework for quantitative 
interpretation of observational data in terms of turbulent properties. 

The observational data has been obtained for the variety of lines, mapping different 
regions of our galaxy (or its satellites) and different physical conditions - 21 cm diffused 
HI in outer parts of the Galactic disk I21I1 . towards the Galactic center 12211 and in 
Small Magellanic Cloud [23], CO lines in molecular clouds l24ll . H a lines in Reynolds 
layer ll25ll . In all the cases the reported index of the power spectrum of turbulence 
is close to « ~ —2.7. This is really remarkable given that these results are obtained 
for quite different observables and that n = —2.7 does not correspond to Kolmogorov 
ll26ll cascade, the one case when some universality may be expected. To complicate the 
matter, some lines are optically thin, given view for the full depth of the emitters along 
the line of sight, while for the other (definitely for CO, but also for HI viewed towards 
the Galactic center) self absorption is important. Here, using the formalism of VCA, we 
will demonstrate how the universal power spectral scaling in the data may be explained. 



STATISTICAL DESCRIPTION OF THE TURBULENCE 

Stochastic density and velocity fields 

The main feature of the stochastic random fields arising as a result of turbulent 
processes is that they are correlated. The central object which describes their properties 
is the correlation function £(r) (or where more suitable - the structure function d(r)). 
Related and often less limited descriptor is the power spectrum P(k), which is a Fourier 
transform of the correlation (structure) functions. 

We shall restrict our discussion to the case of the turbulence that is statistically 
homogeneous and isotropic in three dimensional position space (l7l| . 1 in which case 
correlation functions depend only on separation distance r between two positions in 
space. In particular, for the density field p(x) we define the correlation function 

§(r) = §(r) = (p(»)p(x + r)> . (1) 



In the presence of magnetic field MHD turbulence becomes axisymmetric in the system of reference 
related to the local direction of magnetic field. The Goldreich & Sridhar [7] model of incompressible 
turbulence prescribes the Kolmogorov scaling of mixing motions perpendicular to magnetic field lines 
with D(r) °< r 2 / 3 and a different scaling along the magnetic field. Further research in Ulal has shown 
that the basic features of the Goldreich-Sridhar turbulence carry over for Alfvenic perturbations to 
the compressible regime. Observations, however, are usually unable to identify the local orientation of 
magnetic field and deal with the magnetic field projection integrated over the line of sight. As the result 
the locally defined perpendicular and parallel directions are mixed together in the process of observations 
fToll . There is some residual anisotropy, but this anisotropy is scale-independent and is determined by 
the rate of the meandering of the large scale field. So from the observational point of view the picture of 
the isotropic turbulence remains to some extend applicable. The spectra of intensity fluctuations obtained 
from the Goldreich-Sridhar turbulence and the isotropic turbulence are similar. 



or, alternatively, the structure function 

d(r) = ((p(x + r)-p(x)) 2 ) 



(2) 



Statistical descriptors can often be assumed to have power-law dependence on scale 

r~ r , where 7 can be both positive and negative. The correlation function £ (r) is an 
appropriate choice when stochastic inhomogeneities have more power on small scales, 
which corresponds to 7 > 0, while when the power is concentrated on large scales and 
7 < 0, the structure function should be used. With this substitution in mind both cases 
can be treated similarly in the most interesting regime of small r. 

For power law correlations the power spectrum P(k) = / d N re' kr ^ (r) also has power 
law form P{k) k", where n = 7 — N. It is important to keep in mind that the relation 
between spectral index n and correlation scaling 7 involves the dimensionality of space 
N and is different for three dimensional fields and two dimensional maps. 

To be exact, in the case of density it is more appropriate to assume a power law scaling 
for the density fluctuations around the mean 8p = p — p rather than for the density field 
itself, so that 

S(r)=p 2 [l + (r/r )-r] . (3) 

For example, in Kolomogorov turbulence the medium is incompressible, however the 
passive tracers of the flow, which we may observe as emitters develop passive scalar 
density inhomogeneities which scale as r 2 / 3 , i.e 7= —2/3 in our notations. 

Whether the medium is incompressible or not, the stochastic motions are the essence 
of turbulence. An isotropic random velocity field u(x) is fully described by the structure 
tensor (AwjAw,-), which can be expressed via longitudinal Dll and transverse D^n 
components IU7I1 

(AuiAuj) = iPuir) -D NN {r)) ^ +D NN (r)8 ik , (4) 

where 5^ equals 1 for i = k and zero otherwise. We need only projection of the velocity 
structure function onto the line of sight, which we identify with z direction D z (r) = 
(AuiAuj)ziZj. For the power-law regime 

D z (r)~CV" 2 . (5) 

In Kolmogorov turbulence velocity scales the same way, m = 2/3, as the density of the 
passive scalar. Thus the Kolmogorov spectrum index is — 1 1/3. In turbulence literature 
the energy spectrum E(k) = 4nk 2 P(k) is usually used. In this notation the Kolmogorov 
spectrum is E(k) ~ k~ 5 / 3 , often referred to as —5/3 law. 



There is a residual dependence of D z on the angle cos 9 = r • z which differs for solenoidal and potential 
flows, see fl9l 12011. VCA formalism is shown to be insensitive to it, so an important problem of separating 
solenoidal and potential components of the velocity has to be addressed by combination of techniques, 
possibly by combining VCA and velocity centroids |27|]. 



Statistical description in Position-Position- Velocity space 



One does not observe the gas distribution in the real space galactic coordinates xyz 
where the 3D vector x is defined. Rather, intensity of the emission in a given spectral line 
is defined in Position-Position- Velocity (PPV) cubes towards some direction on the sky 
and at a given line-of-sight velocity v 3 . In the plane parallel approximation the direction 
on the sky is identified with xy plane where the 2D spatial vector X is defined, so that 
the coordinates of PPV cubes available through observations are (X,v). The relation 
between the real space and PPV descriptions is defined by a map (X,z) — > (X, v). 

The central object for our study is a turbulent cloud in PPV coordinates. The number 
of particles per unit volume dXdv is given by the PPV density p s (X, v), which statistical 
properties depend on the density of gas in real galactic coordinates, but also on velocity 
distribution of gas particles. Henceforth, we use the subscript s to distinguish the quan- 
tities in (X, v) coordinates from those in (X,z) coordinates. We shall always assume 2D 
statistical homogeneity and isotropy of p s (X, v) in X-direction over the image of a cloud. 
However, homogeneity along the velocity direction can only be assumed after additional 
considerations, if at all. Naturally, there is no symmetry between v and X. 

In lEoll we discussed in detail the derivation of the correlation function in PPV space 



&(*,vi,v 2 ) = (p,(Ri,vi)p,(R 2 ,v 2 )>, R= |Ri-R 2 | 



(6) 



in a general case of matter concentrated in a finite cloud and possibly subjected to 
coherent non-random flow, for example rotation. The main result is that for small scales 
the complications arising from the finiteness of the matter distribution or its non-random 
motion can be neglected and 



&(R,v) 



dz 



A(r)+2j8]V2 



exp 



2(D z (r) + 2/3) 



V = Vi -v 2 



(7) 



where /3 = kgT /m is related to the thermal velocity of atoms. How small should the sep- 
aration R be ? In case of a cloud of size S, obviously we need R <C S to neglect boundary 
effects. In case of the coherent flows, it is the magnitude of the shear introduced by the 
flow relative to the magnitude of the stochastic motions that is important. Exact crite- 
rion is R <C [f 2 C] 2 "' EH, if we parameterize the line-of-sight variation of the coherent 
velocity as v co h = f~ l z. In case of our Galaxy disk rotation, « 14 km/s/kpc, while 
relative turbulent motions reach 30 km/s at few tens of parsecs separation. In Ill9ll we 
have concluded that for scales less than 100 pc or perhaps even a bit more, the rotation of 
the disk is not important. This argument, on the other side, demonstrates that HI velocity 
is a poor indicator of the real line-of-sight position in PPV cubes. 

The equation © is an important and rather universal result. However, it contains 
several approximations which entered the derivation and have to be kept in mind: 

• We have assumed that the underlying turbulent velocity obeys Gaussian statistics. 
No such assumptions have been made for the density, however. 



3 All velocities are the line-of-sight velocities, we omit any special notation to denote z-component. 



• Perhaps most importantly, we have assumed that the velocity and density are 
uncorrelated. It is definitely true if they are taken at the same point due to vector 
nature of velocity field, but needs not be accurate at finite separations. An example 
of a process which may produce correlations is the self-gravity of the gas, which 
will point some vel ocity vectors towards overdensities. This assumption was tested 
numerically in [28, 29]. Correlations have to be quite high to change our results, 
and that level of correlation has not been observed. 

• At the final stage, we assume homogeneity along the velocity dimension as well 
as in two spatial directions on the sky. For this to be applicable in practice, one 
may have to subtract the average spectral line profile and work with fluctuations of 
intensity. 

• The formalism, leading to the equation Q has to be extended if the temperature of 
the gas varies from point to point and the turbulence is near subsonic. 

The correlation in PPV space reflects, as expected, both underlying inhomogeneities 
of the matter through t, (r) and its velocity through D z (r). Important feature of the PPV 
space is that p i5 (X,v) exhibits fluctuations even if the flow is incompressible and no 
density fluctuations are present. Indeed, when one substitutes the expanded expression 
©, E, (r) = p 2 + p 2 (/o/r) r , into eq Q, both terms give rise to non-trivial contributions 
toij/(R,v). 



FROM THE LINE INTENSITY FLUCTUATIONS IN VELOCITY 
CHANNELS TO THE CHARACTERISTICS OF A TURBULENT 

CASCADE 

Intensity 7 V (X) in the spectral line measured at the velocity v in the direction X is 
obtained by integrating the standard equation of radiative transfer (Spitzer 1978) 

dl v = -g v I v dz + j v dz (8) 

along the line of sight. In the case of self-adsorbing emission in spectral lines the 
coefficients are proportional to first power of density: 

g v (z) = a(z)p(z)0 v (z) , 

j v (z) = ep{z)Uz) , (9) 

where (j) v (z) describes the velocity distribution of atoms at the position z along the line 
of sight, and all quantities have implicit dependence on the sky direction X. Obvious 
relation / dzp(z)(j) v (z) = p s (X,v) links us to the formalism of PPV density. 

A solution of the radiative transfer equation if no external illumination is present and 
absorption coefficient a can be taken as constant (the latter is the essence of the Sobolev 
approximation) can be written in a compact form 12011 

(10) 



«X) = - 



1 _ e -ap*(X,v) 



In the case of vanishing absorption, the intensity is given by the linear term in the 
expansion of the exponent in eq. ([TOl 

/ v (X)=ep,(X,v) (11) 

and reflects the PPV density of the emitters. If, however, the absorption is strong, 
the intensity of the emission is saturated at the value e/a wherever p s (X,v) ^> l/a. 
Identification of the low contrast residual fluctuations may be difficult in practice. 
In observations we register the intensity /c(X), integrated over velocity channel 

/CO 
dvW(v-v*)I v (X) (12) 
-co 

which width and shape is described by the window function W(v). The properties of the 
window function are restricted, first of all, by the experimental setup. For instance for 
CO lines studies [30, 24] integration is typically being performed over the whole line, 
W(y) = 1. On the other hand, measurements in 21 cm HI line are performed in velocity 
slices of PPV data cube (channel maps), which corresponds to W(y — v*) strongly peaked 
at a particular velocity v*. 

When velocity channels are narrow we have an opportunity to synthetically increase 
the width of the channels by combining the data from the adjacent ones. Within the 
VCA technique we point out that by varying the width of synthetic velocity channels 
one can separate statistics of turbulent velocities and density inhomogeneities of emit- 
ting medium. The minimal width of the velocity channel is determined by the resolution 
of an instrument. However, thermal motions of the gas lead to the smearing of intensity 
fluctuations similar in effect to finite resolution. Qualitatively, one can think of the min- 
imal width of the channels available for VCA as given by a convolution of instrumental 
and thermal effects. 

Main statistical quantity that we measure from two dimensional intensity maps /(X) 
is the structure (correlation) function 

0(R) = ([/c(Xi)-/c(X 2 )] 2 ), R = Xi-X 2 , (13) 
S(R) = (/ c (Xi)/ c (X 2 )), (14) 

(or, alternatively, 2D power spectrum). Here I shall proceed to describe how underlying 
three dimensional statistics are reflected in this function. 

Optically thin lines 

From the point of view of VCA approach, the most informative data comes from 
optically thin lines, such as, for example HI 21 cm. In the limit when the absorption can 
be neglected line intensity contains direct information about PPV density of emitters 
(fTTTi and the correlation function is simply 



0(R) = e 2 / dviW(vi) / dv 2 W(v 2 )[d s (R,v h v 2 )-d s (0,v h v 2 )} , (15) 



S(R) = e 2 JdviW(v 1 )Jdv 2 W(v2)U^n,V2). (16) 

Let us assume that after subtracting the mean line profile, the fluctuations of PPV 
density are statistically homogeneous not only in X, but in velocity direction as well 4 . We 
shall write out in detail the correlation function as leading to more compact expressions, 
but present the final results in terms of more appropriate @(R). From the eq. 0, in the 
expanded form, 



Z(R) oc e 2 f dz^(r) D z ( r y 1/2 f dv W 2 (v) exp 

J — CO J 



V 2 



2D z (r) 



(17) 



where v is the difference in velocities between two emitters, v = vi — v%, v + = vi + V2 and 
the effective channel window W 2 (v) is a square of the experimental window convolved 
with the gaussian that describes the thermal velocities 



W 2 (v) = -J= J dye V J dv+ Wc(v+-y)W c (v + +y) . 

For convenience let us parameterize the effective window with a gaussian shape with 
the width AV, W 2 = AV~ l exp[— v 2 / (2AV 2 )}. The exact shape, of course, depends on 
the properties of the instrument. 

Equation ([T71) provides the main theoretical foundation for VCA of optically thin 
lines. It describes two fundamental regimes: 

thick channel The effective channel width AV is larger than the turbulent velocities at 
the scale R, AV 3> D Z (R) 1 / 2 . In this case we collect inside our channel essentially 
all emitters, however scattered along the velocity axis they are by the turbulence. 
Essentially, W 2 ~ 1 and 

/•OO 

Z(R)oc / dz^(r) . (18) 



Information about the turbulent velocities is erased, but one can recover information 
about density inhomogeneities. 

thin channel , AV D z (R) l l 2 . Turbulent velocity differences are important, many 
emitter pairs close along the line of sight acquire disparate velocities and do not 
contribute to the correlation of sources restricted to the narrow velocity channel. 
The amplitude of correlation decreases, but the most important outcome is that this 
decrease is scale dependent and we have a change in scaling law which reflects the 
underlying velocity statistics. Indeed, in this case 

/oo 
dzS(r)D z (r)- 1 / 2 . (19) 
-oo 



For an infinite emitting medium homogeneous turbulence produces a homogeneous image in the velocity 
space. For a finite emitting cloud the PPV image is approximately homogeneous over velocity separations 
much less than the Doppler line-width. If we use the structure function, we do not need to explicitly worry 
about subtracting the mean line profile, since its integrand does not depend on it. 



TABLE 1. Scaling of emission intensity in optically thin lines 



%{R) 



P V {K) 



Pp{K) 



Thin channel, < CR m 
Thick channel, > CR m 



OC 








oc 



jr— 3+m/2 




Interpretation of thick channel data contains no surprises, reflecting projected density 
inhomogeneities of the gas. However if measurements are done in a sufficiently thin 
channel (and the criterion depends on scale between two lines of sight), the fluctuations 
of intensity may come from two sources. Recall that £ (r) °c 1 + (r/ro)~ r , therefore 



Z(R)oc dzD z (r)- l ' 2 + dz{r/ro)-7 D z (r)- l ' 2 °c 1 +Z V (R) +Z p (R) . (20) 



The first term describes the intensity fluctuations arising just from turbulent motions of 
emitters 5 . The second is the velocity modified scaling of the underlying density fluc- 
tuations. In Table [T] we summarize the predicted scaling for the optically thin emis- 
sion intensity in the channel maps, giving both the structure functions, and the the 
correspondent two dimensional power spectra P V (K) / JR E V (R) e' KR and P p (K) 



Thus, we are in a position when velocity information can be obtained from thin chan- 
nel data, while increase of channel thickness will determine the density distribution. 
If the density scaling is steep, n < —3, as is the case of emitters, behaving as passive 
scalars in incompressible Kolmogorov cascade, thin channel measurements are veloc- 
ity dominated and exhibit spectra oc ^- 3+m/2 ^ K 21 if the stochastic velocities are 
Kolmogorov, m = 2/3. 

It is a question whether one can obtain thin channels in a particular observations. The 
criterion is scale dependent, for sufficiently small separation between the lines of sight, 
the effective width of the channel is always large, but the needed larger separations may 
be more difficult to measure. Thermal effects also increase AV. Thus, observing in thin 
channels may require special design and be achievable only for cold gas. However, as we 
have shown in lfl9ll . the HI observations, in particular lEHl23ll . are, effectively, the thin 
channel measurements. Thus we may interpret the observed slope of —2.7 as a signature 
of the Kolmogorov stochastic velocities, but we cannot exclude an alternative that the 
underlying density inhomogeneities have enhanced small scale power and n « — 3. 



Relation to the underlying turbulence of the intensity in spectral lines that exhibit self- 
absorption is much more involved 12(11 . If the absorption is strong, the intensity of the 



5 this term contains constant which corresponds just to mean intensity in the map. It has to be regularized 
out to make sense of the integral, which is done automatically if structure functions are used. 
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Effects of self-absorption on spectral line statistics 




0.001 



0.010 



0.100 



R/S 



FIGURE 1. Structure function of intensity fluctuations from Kolmogorov turbulent motion in case of 
weak self-absorption in the spectral line. 

emission dTOl) is saturated at the value e/a wherever p 5 (X, v) ^> 1/a. Identification of 
the low contrast residual fluctuations may be difficult in practice. 

However the study of such lines is crucial, in particular for understanding molecular 
clouds. To make a progress, we look for the regimes when intensity spectra can still 
exhibit scaling response to the power-law turbulent cascade. In H20I1 we have found that 
in the case of weak but not zero absorption, the universal scaling solution Q>{R) R (this 
corresponds to K~ 3 for the power spectrum) arise over the range of intermediate scales. 

Let us return to the structure function S>{R). One can always consider sufficiently 
small scales such a 2 [d s (R,v) — rf 5 (0,v)] <C 1. If this holds, we argue in [20] that 



provides a good (although not fully rigorous) approximation to intensity structure func- 
tion. Although this expression looks similar to eq (IT6l) for the optically thin case, the 
conditions for its validity are more relaxed than a limit a — > 0. 

The most important effect induced by absorption is retained. It is an additional expo- 
nential down- weighting of the contribution from the points with large velocity separation 
v in a manner which itself depends on the turbulence statistics. This has an effect of the 
stochastic velocities having an impact even on the fully integrated over the frequency 
lines. 

Here we limit ourselves to the short discussion of this effect, considering densit y in - 
homogeneities to be subdominant in PPV statistics. For more general issues see 112011 . 
Numerical calculations of the S)(R) for Kolmogorov underlying velocity given in Fig- 
ure [T] show that at large scales R the effects of absorption become highly nonlinear 
(indeed, the approximation d2TT) breaks down), correlations are erased and the structure 
function flattens. At shorter scales, where one may perhaps expect the regime of ab- 
sorption behaving like like a thin velocity channel, the structure function exhibits the 




(21) 



TABLE 2. Scalings of structure functions of intensity fluctuations arising from velocity fluctu- 
ations for the power-law underlying 3D velocity statistics. In the strong absorption regime @(R) 
does not follow a simple power-law. The regime of thin channel is not realized for m > 2/3, 
including the Kolmogorov case. 

Scale Range Intensity Scaling Regime 

R/S < (fi/D z (S)) l / m velocity effects erased (subsonic regime) 

/ 2 \ '/'" 

R / S < TCTsT @( R ) x R l+m/2 (transition to thick slice) 



DAS), 
-. \ 1/m 



R / S< [d^)J 2>{R)^R x (intermediate scaling) 

2 \ !/'" / 2 \ 2 /( 2 ~ m ) 

■jgfa) <R/S<(rffsj) &{R)ocR 1 -"'/ 2 (thinshce) 

- v2 s 2/(2-m) 

J < R/S not a power law (strong absorption regime) 



R 1 slope instead, which we found to be independent on underlaying velocity scaling 6 . 
At even shorter scales the absorption becomes negligible and the ordinary thick channel 
regime takes over. 

What determines the scales of transition ? We may think of absorption as providing a 
window with the width Av ab defined by 

a 2 d s (0 1 Av ab ) = \ . (22) 

Only the matter with relative velocities less than Av ab contribute to observed intensity in 
a correlated fashion. Smaller is the Av ab , stronger are the absorption effects. This critical 
relative velocity depends on the mean density in PPV space, the absorption coefficient, 
and the slope of the PPV structure function. How important the absorption is for intensity 
correlations at a scale R depends on comparison between Av ab and rms turbulent velocity 
at this scale, D Z (S) l ' 2 (R/ 'S) m l 2 where S is the size of our cloud. If Av a j, is larger than rms 
velocity differences between emitters, its effect is gradually diminished. This translates 
into absorption being less important on small scales and more important on large scales. 
A detailed account of different regimes and transition scales is given in Table El 

In view of this analysis CO measurements of 12411 which show shallower, —2.7, slope 
instead of —3 may be taken to indicate that in molecular clouds the density inhomo- 
geneities have enhanced small scale power with n > — 3 and dominate the velocity con- 
tribution. However to make definitive statement one needs to make sure that the differ- 
ence of 0.3 in the power spectrum can not arise from the noise in the data. 8n = 0.3 or 
better seems to be an accuracy we need to have to distinguish between different theoret- 
ical possibilities discussed in this paper. 



The asymptotic analysis in 1 20] shows how this behavior arises 



SUMMARY 



We have reviewed the formalism of Velocity Channel Analysis for optically thin lines 
and its extension to the lines with self-absorption. We demonstrate that by observing 
optically thin lines from cold gas in sufficiently narrow (thin) velocity channels one 
may recover the scaling of the stochastic velocities from turbulent cascade, in particular, 
Kolmogorov velocities give K contribution to the intensity power spectrum. This 
matches the observational data in HI, both in our Galaxy and in Small Magellanic cloud. 
Synthetically increasing the channel thickness separates out the underlying density 
inhomogeneities of the gas. An attempt to apply this technique to HI lines in SMC [12311 
indeed showed the expected steepening of the spectrum with the channel thickness. 

Effects of self absorption, on the other hand, retain the velocity signature even for 
integrated lines. As a result, intensity fluctuations from velocity tend to show universal 
scaling of the power K ^ over the range of scales. Observed shallower spectra n m 
—2.6 in CO lines therefore may indicate that real density spectrum is also shallow 
n ps —2.6 and dominates the signal. 

Progress in interpretation of the data require the knowledge of the power spectrum 
index better than with an accuracy of 0.3 and developing of the accurate theoretical and 
numerical models to match such accuracy in different regimes. 
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